Role of anisotropy in the spin-dimer compound BaCuSi2 06 
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We present results of magnetisation and electron paramagnetic resonance experiments on the 
spin-dimer system BaCuSi2 06. Evidence indicates that the origin of anisotropic terms in the spin 
Hamiltonian is from magnetic dipolar interactions. Axial symmetry-breaking is on a very small 
energy scale of < 11 mK, confirming Bose Einstein condensation critical scaling over an extended 
temperature range in the vicinity of the quantum critical point. 



A field-tuned quantum critical point (QCP) separates 
the low field quantum paramagnetic phase in spin-dimer 
compounds from the magnetically ordered phase in high 
magnetic fields. The XY antiferromagnetic ordered state 
of these compounds corresponds to a Bose-Einstein con- 
densate (BEG) in the absence of axial symmetry breaking. 
To ascertain the extent of the BEC universal region in 
each spin dimer material therefore, it is vital to iden- 
tify the nature and size of any anisotropic terms that 
may be present in the Hamiltonian In order to ob- 
serve BEC critical scaling, any anisotropic terms must 
be small enough that the universal region is at least a 
decade in reduced field. Anisotropic terms may be di- 
rectly measured by means of Electron Paramagnetic Res- 
onance (EPR) , which has revealed for example in the pro- 
totypical spin dimer cornpound TlCuCls (triplon band- 
width ~i mev) mmiiii, staggered Dzyaloshinskii- 

Moriya (DM) terms on the order of 0.1 meV j^. In this 
report, we discuss the size and role of anisotropy in the 
spin dimer compound BaCuSi206 1^. Hol llH . 

A first approximation to the equivalent magnetic lat- 
tice corresponding to the room temperature tetragonal 
body-centered BaCuSi206 lattice 0| is the isotropic spin 
Hamiltonian: 
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An antiferromagnetic (AF) exchange constant J > 
couples each pair of Cu^+ ^ ~ ^ spins within ver- 
tical spin dimers on Cu2Si40i2 layers (where i is the 
index of the dimer, and Sn and 8^2 are the two spins 
that form the pair i.) Intra-layer (J') and inter-layer 
(J/) AF exchange constants couple nearest neighbour- 
ing spin dimers on the lattice, indexed by Cq, = {x,y}, 
= {zdzx/2±?//2} as defined in the high-temperature 
structure [l2|. For an s = i system, energy levels cor- 
responding to the isotropic Hamiltonian are a ground- 
state singlet, and three degenerate triplet excited states, 



separated from the groundstate by the spin gap A. An 
applied magnetic field H introduces a Zeeman splitting 
term in the Hamiltonian, scaled by the ij-factor gg, cor- 
responding to the orientation of H. Experimental results 
reveal a critical magnetic field i?ci ~ 23.5 T at which the 
system orders for i? || c, and comparison of the measured 
phase boundary with the model Hamiltonian in Eqn. ^ 
gives values of J = 4.45 meV, J' = 0.51 meV, and an 
even smaller value of J/ 0, 0| . 

The above description neglects anisotropic exchange 
interactions which can reduce the symmetry of the or- 
dered state. In the more general spin Hamiltonian, the 
bilinear spin term • Sj = Si I Sj in Eqn. ^ is replaced 

by the exchange interaction j=i ^i^jj^j; where 



^as TgjYi, (2) 
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Ts is a scalar, while the antisymmetric term Tas of the 
form D • [si X Sj], and the traceless symmetric term 
of the form s^Fsj mix spin components 0|. The DM 
interaction arising from the spin-orbit coupling is anti- 
symmetric in nature, leading to singlet-triplet {ST) mix- 
ing to the lowest order of perturbation. Magnetic dipolar 
interactions, however, belong to the traceless symmetric 
category. Terms of this nature lead to intra-triplet ma- 
trix elements entering to the lowest order of perturbation, 
but ST mixing to higher orders of perturbation only for 
H oriented away from the dipole vector. An additional 
possible source of anisotropic exchange interactions is the 
reduced symmetry of the BaCuSi206 lattice at low tem- 
peratures (recent structural analysis has revealed a weak 
orthorhombic structural transition at ~100 K, accompa- 
nied by an incommensurate lattice modulation |l5j|l. 

An applied magnetic field with a component perpen- 
dicular to cither the DM or dipolar vector would result 
in spin non-conservation due to mixing of singlet and 
triplet energy levels with inequivalcnt Sz- The symme- 
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try of the Hamiltonian would therefore reduce from rota- 
tional invariance (U(l)) to the discrete Z2 group, result- 
ing in a magnon spectrum in the ordered phase with a 
finite gap to the minimum, rather than a gapless Gold- 
stone mode. The extent of symmetry breaking, however 
depends on the order of perturbation at which spin non- 
conserving terms enter the spin Hamiltonian; it is there- 
fore important to identify the origin of any anisotropic 
interaction in the Hamiltonian. Previous measurements 
of a BEG critical exponent in BaCuSi206 down to 0.03 
provide empirical evidence for the absence of 
U(l) symmetry-breaking terms down to this energy scale 
for H\\c. In this rapid communication, we discuss more 
direct experimental evidence that estimates the size of 
any axial symmetry-breaking due to anisotropic exchange 
terms in the spin Hamiltonian. 

Single crystals were grown using the slow cooling flux 
technique. Polycrystalline BaGuSi2 06 precursor was 
synthesized by a solid state reaction of BaGOa, Si02, 
and CuO between temperatures of 900 and 1050°G in 
flowing oxygen, with repeated regrinding. Single crystals 
up to 1 g were grown by heating a 2:1 molar mixture of 
ground polycrystalline material and LiB02 flux in a plat- 
inum crucible to 1000°G in air, followed by slow coohng 
to 875° and decanting by centrifuge. 

EPR measurements were performed for a range of fre- 
quencies between 26 GHz and 660 GHz for both H \\ c 
and H J- c. At elevated temperatures, population of the 
Zeeman split triplet states results in intra-triplet transi- 
tions, observed as sharp dips in the transmission through 
a microwave cavity (Fig. ^) . At the highest tempera- 
tures, a single EPR peak is observed, corresponding to 
the degenerate transitions between triplet levels (in the 
absence of anisotropic terms). The intensity (integrated 
area under the peak) diminishes rapidly upon reducing 
the temperature below 30 K (Fig. ^), as the triplet 
states depopulate. A residual intensity can be observed 
to persist, and even grow, at the lowest temperatures, ex- 
hibiting hyperfine splitting (into 4 distinct peaks). The 
4-peak pattern is consistent with hyperfine coupling as- 
sociated with a localized s — and a lone Gu nucleus 
(/ — |). From fits to the temperature dependence of 
the total integrated intensity (Fig.^), we conclude that 
the central portion of the spectrum observed below about 
5 K is dominated by a small concentration (0.5(1)%) of 
isolated paramagnetic Cu^"*" defects (s = 5) in agree- 
ment with estimates from DG magnetic susceptibility 
(x — rUz/H) measurements of the same sample made 
using a SQUID magnetometer (Fig. ^). Since the hy- 
perfine splitting has the same g factor and anisotropy as 
the lattice, the associated Gu^+ defects likely correspond 
to singly-occupied dimers. Perhaps the most important 
feature in the data is the appearance of distinct, albeit 
weak shoulders which split from the central peak at low 
temperatures (dashed lines in Fig. ^). The vanishing 
intensities of these peaks as T — > indicate that they 
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FIG. 1: [Colour online] (a) EPR spectrum measured at dif- 
ferent T for H\\c and H c using a resonant frequency of 
51.8 GHz. For clarity, data for different T have been offset. 
Weak shoulders appear and split from the central EPR peak 
upon reducing the temperature (see dashed lines and main 
text for explanation). The central peak also exhibits hyper- 
fine splitting for -ff || c (see main text). The inset shows a 
schematic of the possible transitions between different spin 
states, (b) temperature dependence of the magnetic suscepti- 
bility for a field of 5000 Oe aligned parallel and perpendicular 
to the crystalline c-axis. (c) measured temperature depen- 
dence of the integrated EPR intensity for H\\c. Solid lines in 
panels (b) and (c) show fit to isolated dimer model including 
a Curie term to account for paramagnetic impurities. 

involve EPR transitions within the triplet state. 

The frequency and field orientation dependences of the 
split shoulders are displayed in Fig. [3 The splitting sug- 
gests an anisotropic zero-field-splitting (zfs) interaction 
in the spin Hamiltonian, as confirmed from the frequency 
dependent measurements for H \\ c, i.e the splitting {5) 
is field-independent. The field orientation dependence of 
the splitting can be fit by the equation: 

(5= ^i?(3cos2 - l)|^lBge (5) 

where 9 is the angle between H and the vector between 
the intra-dimer Gu^+ sites (i.e. the c-axis). One possi- 
ble source of such an angular dependence is a symmet- 
ric anisotropic interaction (Eqn. ^ in the Hamiltonian 
of the form DS^, where S is the total spin (= 1) of 
the dimer and D — T^z- The value of the anisotropic 
parameter D is found from the angle dependence to 
be 0.091 (3) K. Similar shoulders have been observed 
in Ref. [13, and the angle dependence of the shoul- 
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ders fit to D K. A more direct measure of the 

anisotropic parameter involves a determination of the 
zfs of triplet levels by extrapolating the frequency de- 
pendence to zero field (Fig. |2J). By this method, zfs of 
1.8(3) GHz [0.09(1) K] and 2.1(1) GHz [0.10(1) K] are 
obtained, in agreement with the value from the angle de- 
pendence. We note also that fits to the central high tem- 
perature peak yield extremum values for the g-factor of 
g\\ — 2.307(3) and g±^ = 2.057(3), in excellent agreement 
with susceptibility and X-band EPR data • 

The observed zfs is consistent with either an antisym- 
metric DM interaction (Eqn. OJ, with D || c J^] and 
D| = 4.8(2) K or a symmetric dipolar intra-dimer in- 
teraction. A DM interaction would mix singlet {S) and 
triplet ^2 = (To) states for i/||c to lowest order - which 
although not U(l) symmetry-breaking, would allow a 
significant EPR transition between the S and triplet 
= -1 (r_) states (labelled ST in Fig. [TJ. Simula- 
tions which include an antisymmetric anisotropy predict 
similar intensities for the triplet transition (labelled Rl 
and R2 in Fig.^ and the high- field ST transition at 6 K. 
In contrast, for H\\c, an ST EPR transition is forbidden 
by symmetry for an intra-dimcr dipolar interaction, and 
would be negligibly small due to higher order mixing of 
the S and T_ states for an inter-dimer dipolar interac- 
tion. For H J- c, the extent of mixing of the S and 
T_ states (which would break U{1) symmetry) would be 
considerable for a first order DM process, but very weak 
for higher order intra- and inter-dimer dipolar processes. 
Therefore, the DM interaction would lead to significant 
non-linear field-dependence of the intra-triplet transition 
frequency for fields approaching Hd, whereas the dipo- 
lar interaction would not. High- field EPR measurements 
reveal no sign of ST transitions in the vicinity of Hd, or 
any unusual behaviour of the intra-triplet EPR transi- 
tions for ± c up to 21.7 T (Fig.lSJ, strongly indicating 
that the observed low-temperature EPR splitting arises 
from dipolar interactions which are symmetric in nature. 
Furthermore, as we describe below, both the magnitude 
of the zero field splitting and its T-dependence are con- 
sistent with dipolar, and not DM interactions. 

The dilute concentration of well separated thermally 
activated triplets at low T, H provides the ideal con- 
ditions to observe the effects of magnetic dipolar inter- 
actions in BaCuSi206 . The magnetic lattice comprises 
closely spaced pairs of Cu + ions (r|| = 2.74 A intra- 
dimer spacing at room temperature) which are well sep- 
arated from each other (r_L = 7 A between dimers). 
The zero-field dipolar splitting of triplets is given by 
Y^(25^ + gD^i^g, which has the value 0.113 K for 

r|| = 2.74 A. This is remarkably close to the measured 
zfs ~ 0.1 K, suggesting that the origin is indeed dipolar 
interactions. The reason for the collapse of the dipolar 
splitting at higher temperatures can be understood in 
terms of 'exchange narrowing'. At low T, the triplets 



120 




0.0 0.5 1.0 1.5 2.0 2.5 

FIG. 2: [Colour online] Magnetic field dependence of experi- 
mentally measured intra-triplet splitting (labelled Rl and R2 
following Fig. for H\\c. The solid lines are linear fits to 
the intra-triplet splitting as a function of magnetic field, ex- 
trapolated down to zero magnetic field. The inset shows the 
angular dependence of this splitting measured using a reso- 
nant frequency of 51.8 GHz at a temperature of 6 K. The solid 
lines show a fit to Eqn. |^ with D = 0.091 K. 

700 1 -] 

^ H inclined 64° 

3-1 Ann 

O away from c 

500 - 

1 "^^^ ^ R1&R2 

^ K '^=.^^ ^ 

12 14 16 18 20 22 

FIG. 3: High field linear field dependence of the intra-triplet 
Zeeman splitting with H inclined with respect to [001]. 



are dilute and long lived on EPR time scales (^ I/./)- 
The zfs results from the anisotropic dipolar field that 
each spin within a dimer experiences due to its pair 
[oc 3(si • z){s2 ■ z) — {si ■ 82)]. However, as T is raised, and 
more triplets are excited, any given Cu spin will expe- 
rience strong fluctuations in the local dipolar fields due 
to the exchange-induced co-fiipping between neighboring 
dimers of opposing spin projection. As more triplets are 
excited, such co-flipping leads to faster fluctuations of the 
local dipolar fields until they are eventually averaged out 
on EPR time scales, and the dipolar splitting vanishes. In 
fact, most of the linewidth observed in these experiments 
can be attributed to nuclear and dipolar spin-spin inter- 
actions (both intra- and inter-dimer). As T is raised, 
exchange averaging leads to a gradual reduction in the 
second moment of the dipolar field distribution, and to a 
narrowing of the spectrum. We do indeed observe a fur- 
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ther narrowing of the spectrum at higher T (not shown) , 
which has also been observed in Ref . |3 . 
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FIG. 4: [Colour online] Magnetisation measured as a func- 
tion of H up to 40 T. The applied field H\\c and H _L c 
(sr|l = 2.31, gx_ = 2.05) is represented as an equivalent energy 
E = iJ.ofJ.BgeH meV on the lower a;-£ixis, such that both 
orientations share the same axis. The upper a;-axis shows H 
along the c-axis. 

We also discuss the directional anisotropy in the value 
of g, which arises from crystal electric field splitting of 
the Cu^"*" energy-levels in a tetragonal environment, and 
spin orbit interactions. The anisotropy in magnetic sus- 
ceptibility for _ff ± c and H\\c (shown in Figure re- 
fiects an anisotropy in (7- values. Fits to an isolated dimer 
model using the Bleaney-Bowers equation |1(tI | yield a 
value of J = 4.40(2) meV and values of g± = 2.03(5) 
and g|| = 2.31(5) in agreement with EPR results, where 
II and _L refer to H oriented parallel and perpendicular to 
the c-axis. In addition to scaling the low field suscepti- 
bility, the effect of this g-anisotropy is to scale the value 
of Hci with gg. 

High field magnetisation measurements were per- 
formed in pulsed magnetic fields of 500 ms pulse duration 
up to 40 T. Data were obtained using a wire- wound sam- 
ple extraction magnetometer in which the sample is in- 
serted or removed from the detection coils in situ. Figure 
^ shows the uniform magnetisation (m^) of BaCuSi2 06 
as a function of field for H\\c and H L c. The value of 
ruz for both directions of H has been normalised by 2/ gg 
at 40 K so that the saturation value is 1 hb/Gvl. The 
upturn in niz measured at 0.5 K occurs at Hd, at which 
magnetic ordering takes place. Linear extrapolation of 
niz is used to determine the value of A = ggfisHd for 
H\\c and H L c. The value of spin gap thus extracted is 
A = 3.2(1) meV for H\\c and H J. c. As expected, the 
values of Hd and scale with gg for H along different 
crystal axes. The effect of the anisotropy in g appears 
only in the Zeeman term in Eqn^ indicating that the 
g-tensor is diagonal, and does not lead to mixing of the 



spin components. 

Having carefully examined any possible source of 
anisotropics in the Hamiltonian describing BaCuSi2 06, 
we conclude that the dominant effect arises from symmet- 
ric dipolar interactions, which lead to U(l) symmetry- 
breaking only to second order. The magnitude of U(l) 
symmetry-breaking due to the intra-dimer dipolar inter- 
action is of the order of ~ |-D|J'/J = 11 mK only for 

_L c, whilst of the order of ~ {^)^\D\J'/J = 0.7 mK 
for all H due to inter-dimer dipolar interaction. Indeed, 
the family of spin-dimer compounds is distinct from other 
XY antiferromagnets due to the small size of J' / J and 
f\\/i"±, which 'protect' the U(l) symmetry down to low 
energy scales in the case of symmetric anisotropics. EPR 
measurements on BaCuSi206 confirm that the correspon- 
dence drawn with the BEC universality class is justified 
in a significant region of criticality near the QCP. 
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